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Abstract: We report a study of optomechanical entanglement under 

the drive of one or a series of laser pulses with arbitrary detuning and 
different pulse shapes. Because of the non-existence of system steady 
state under pulsed driving field, we adopt a different approach from the 
standard treatment to optomechanical entanglement. The situation of the 
entanglement evolution in high temperature is also discussed. 
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1. Introduction 

Due to their various interesting properties, optomechanical systems (OMS) are under extensive 
researches over recent years nun . They are regarded as good platforms for realizing high pre¬ 
cision detection 021, e.g. detection of gravitational wave 181191 , ultra-low temperature cooling, 
e.g. cooling the nano-mechanical oscillator to its ground-state CMIl, and macroscopic quan¬ 
tum states, e.g. “cat state” or macroscopic quantum entanglement lT4l422ll . and others. An OMS 
is nonlinear by nature, so it is difficult to find its exact dynamical evolution in the regime of 
strong optomechanical coupling. In the weak coupling regime, the fluctuation expansion around 
the classical steady states is a standard method for studying any feature of an OMS, including 
optomechanical entanglement nma. With the replacements a ^ a + 5d and b ^ b + 5b, the 



original cavity (mechanical) mode is expanded into the sum of the average value d (b), as the 
steady state found from classical equations of motion, and its fluctuation 5d (5b). Then the 
initial Hamiltonian will be linearized so that the quantum Langevin equations about the fluctu¬ 
ations can be solved. The practical use of the fluctuation expansion approach, however, relies 
on the solutions of classical nonlinear dynamical equations, which often take the approximate 
forms. In the cases of continuous wave (CW) drives, such solutions are the classical steady 
states determined by Routh-Hurwitz criterion 

The explicit situations when there is no steady state include the pulse and blue detuned 
CW driving fields. A possible way to deal with these situations is the evolution decomposition 
method 12411251 . which works in any regime without the restriction of the steady state condi¬ 
tion. A number of new phenomena were discovered based on this approach. For instance, the 
entanglement under blue detuned drive is predicted to be higher and more robust against high 
temperature than that under red detuned drive. Moreover, it is found that quantum noise effect 
can be enhanced by drive intensity, and would destroy the entanglement accordingly, resulting 
in sudden death and sudden revival of OMS entanglement II24II25L 

Here we will apply this method to investigate the OMS entanglement created by pulsed 
drives. As it is known to all, a blue detuned drive would “heat” an OMS. A blue detuned CW 
drive, with its longer interacting time with the system, will increase the system temperature 
continually. The decoherence accompanying the rising temperature would diminish and even 
destroy the entanglement from optomechanical coupling. If a pulsed drive is used instead, the 
heating of the system would not be accumulated to significant temperature, and thus have the 
corresponding decoherence safely neglected for the blue detuned drive. This is one motivation 
to consider optomechanical entanglement under pulsed drive. 

OMSs under pulsed drives have been considered for quantum state preparation, cooling, 
quantum entanglement creation, as demonstrated by the recent experiments ||2614T5]| . In 2011, 
Hofer, et al. showed that the entanglement (measured by EPR variance) created by pulsed drive 
of blue detuned (with the detuning Aq = —(Om) can be larger than the upper bound of entan¬ 
glement due to red detuned CW drive ll29ll . Later, He, Kiesewetter, et al. demonstrated the 
existence of quantum steering in pulsed OMS and also predicted that robust EPR steering 
(without requiring a low temperature reservoir) could be achieved by a square pulse iTSll . 

In this paper, we will present a detailed study of entanglement evolution under pulsed drives 
of different detunings and different shapes. With these evolutions, we will show that the created 
entanglement can even last for longer time than pulses themselves, and such entanglement 
is also robust against relatively high temperature when driven by blue detuned laser pulses. 
In addition, if a series of pulses continuously drive an OMS, nonzero entanglement could be 
preserved even when the gap time between the laser pulses is considerable. 

2. System Hamiltonian 

Here we consider an OMS driven by a single laser pulse or a series of laser pulses, which are 
given in the form Y.jE'jit — j ■ with j being the number of pulses, to being the gap 

time between the pulses, and (Oo being the central frequency of the pulsed drive. The function 
E{t) represents the profile of the pulse. Under the pulsed drive, the system Hamiltonian in the 
interaction picture with respect to Hq = (Oca^d + (Omb^b {h= 1 ) is 


Hs(t) = -g{be -f Pe'‘^‘}d^d + i 





j 




( 1 ) 


where g is the optomechanical coupling factor, cOc (co,„) is the cavity (mechanical) oscillation 
frequency, Aq = Wc — CUo is the detuning. In the above equation, the first two terms describe 


the interaction between mechanical mode and cavity mode, while the last terms represents the 
driving on the cavity. 

The system we consider is an open system, i.e. the cavity and mechanical mode will damp 
with the rate K (Ym)- The coupling to the reservoir can be described as the stochastic Hamilto¬ 
nian 1361: 

Hoit) = i{y/Kd'^^cit) + +H.C. (2) 

where i^m) is the stochastic Langevin noise operator. Under these factors, the joint evo¬ 
lution of the system and reservoir manifests in term of the evolution operator t/s(f,0) = 

^g-‘Iodt[Hs{T:)+HD{T:)) ^ Since the coupling and dissipation processes in the evolution are non- 
commutative, it is impossible to get a closed form of this time order exponential. 


3. Factorization of joint system-reservoir evolution operator 

Here we adopt a different approach from the standard fluctuation expansion method to study 
the system evolution. This approach is based on the factorization of evolution operator and is 
valid for arbitrary driving field including the pulsed ones 1241 l25l . The factorizations for the 
evolution operator is given as follows l24l : 




Xe ‘IodT:Hi{T) Left Factorization 

Right Factorization 


where Vj.(f,T) = Texp{—i for k = 1,2. With these factorizations, the evolution 

operator can be simplified, so that it is possible to applied the factorized operators in succession 
to find the system observables that should be obtained by directly acting the original evolution 
operator Us{t,Q) on system operators or system’s initial quantum state. In addition to OMS, 
this method has been applied to other physical systems recently ETHISI. 

The factorization procedure for our concerned system is as follows. At first, we implement 
a right factorization to separate the noise part t/£)(f,0) = Texp{—I /q out of the joint 

unitary operator: 

Us{t,0) =Texp{-/^ dTUDit,T)Hs{T)ul{t,T)} Texp{-/^ c/t//£,(t)}. (4) 


In the main Hamiltonian Uoit, T)Hs{T)U^{t, t), the cavity, mechanical mode operators d, b are 
transformed to 


= e =A(f,T), (5) 

UD{t,T)bUl{t,T) = +limit, t) =B{t,x), 


whsK licit, 1 ) = y/x l[dx'e andn„,(f,T) = are the 

induced quantum noise operators satisfying the following commutation relations: 


rc(f',T) = 

r,„(f',T) = 


[hcit,t'),nlit,x)\ 
[limit,t'),n]„it,x)\ 




(6) 


The transformed system operators therefore satisfy the equal-time commutation relation 
[A(f,, t)] = [Bit,x),&^it,x)\ = 1. After that, we factorize the drive term t/£(T,0) = 


Texp{/o^ dt'Ej{t')A\t to the left side, and the cavity mode in the rest opera¬ 
tor is determined by the transformation 

t4(T,0)A(f,T)t/£(T,0) = A{t,T) + Y, r 

j L 

-f rdtTcit', 

Jo 

= A(f,T)+£D,-(T), (7) 


Finally, we factorize the term Uc{t, t) = Texp{!g lidt' +B\t, A^A(f,f')} 

out of the optomechanical coupling operator from the right side. To the first order of the op¬ 
tomechanical coupling constant g, we obtain the following effective Hamiltonian, 


HoMiT) = + 


^A(f, t)d;. (T) + (f , t)D, (t) + ^ (t)D,(t) 

./' k j k 


( 8 ) 


which is a good approximation of the effective optomechanical coupling Hamiltonian for the 
weak coupling regime of an OMS. The joint unitary operator is therefore factorized as 


Us{t,0) = UE{t,0)UoM{t,0)Uc{t,0)UDit,0). 


(9) 


Since the OMS we consider is an open system, the joint initial state is the following product 
state 


x{0) = p(0) 0 /?(0) = |0), (0| 0 f; ,, |n)„ (n| 0 /?(O), (10) 

n=0 [k+nm) 

where R{Q) is the reservoir state. The operation t/£)(f,0) and UK{t,Q) keep the joint initial state 
X (0) invariant. So the expectation value of a system operator O will be reduced to 

Tr 5 {Op(f)} = Trs{dTr«{t/s(r,0)p(0)®/?(0)t/J(f,0)}} 

= Trs^R{ul^it,O)ulit,O)dUEit,O)UoM{t,O)piO)0R{O)}. (11) 


4. Entanglement calcnlation 

We work in the regime of g <C 1, so an initial Gaussian state will still be Gaussian after OMS 
evolution. For bipartite Gaussian states, the logarithmic negativity of the correlation matrix 
(CM) can well measure their entanglement. The CM is defined as 


y = 


A 



( 12 ) 


with the elements Vij{t) = 0.5(m,M; -f ujUi) — (ui){uj), where u = {xc{t),pc{t),Xm{t),Pm{t)Y■ 
The expectation values in the expression are calculated with Eq. (11). Then the corresponding 
logarithmic negativity is given as ll40lj42ll 


£4^ = max[0, —ln2T7 ], 


( 13 ) 







where 


and 



(14) 


E = detA + detB — 2detC. (15) 

To find the logarithmic negativity experimentally, one should measure the correlations 
for the cavity and mechanical modes. One proposal for doing so is discussed in lfT4l . More 
efforts should be taken to find the feasible ways to realize the measurement of logarithmic 
negativity. 

Since the drive operation t/£(f,0) only displaces the system operators, it will not contribute 
to CM elements and can be neglected in the calculation of logarithmic negativity. Considering 
the expectation value in Eq. (fTTT i. what we should take into account is only the effective Hamil¬ 
tonian in Eq. ([8]l. Under this Hamiltonian the evolution of system operators are determined by 
the following differential equations: 

j 

J 

ge-i’^+r,n){t-r)/2^m„rY^ [D*{T)a+DjiT)a^] 
i 

j 

ge“^'nr-^^i<-^)YJl^D){x)Dk[x). (16) 

j k 

The final term in the second equation can be neglected, since it only brings displacement which 
will not contribute to logarithmic negativity. The terms containing a and b on the right side of 
the equations indicate the beamsplitter (BS) effect between the cavity and mechanical modes, 
while those containing and indicate the squeezing effect which is the main cause of 
entanglement. Both two effects can be enhanced by the higher drive intensity, because each 
terms including the component D{x). Moreover, one sees from the second term on the right 
side of each equation that a higher drive intensity will also enhance the noise effect. 

We rewrite Eq. (fTfil l in the form 


,da 
' dx 


,db 
^ dx 


+ 


-f 

-f 



M{t,x)u, 


where 


M{t,x) 


/ 0 0 
0 0 
—.^sin(cu,„T) —^sin{(OmX) 
^ ^COs{(0mX) ^ cos {(0,nX) 


— J^COs(CUmT) —J^sin(a)mT) \ 

b?COs{(0mX) b?sm{(0n,x) 

0 0 

0 0 / 


(17) 


where ^ = 3ie[P{t, t)] is the real part and = 3e[P(f, t)] is the imaginary part with P{t,x) 






2ge Then the solution of above equation is found as 

i7(f) = g/oS(0 ) 

= (cosh(v4^))/+ ^/'^^' ^plK{t,0))U{0), (18) 

\/'w(f,0) 

where K(t,0) = jQdTM{t,x), and the function m{t,0) is from the relation K^{t,0) = m(t,0)I. 

5. Entanglement under different laser pulses 

5.7. Entanglement evolution under a single laser pulse 

With the solution in Eq. (fTsT l. one can obtain the entanglement evolution under arbitrary 
pulsed drive. The first one we discuss is the Gaussian pulse, whose profile is E{t — jo) = 
jt-j^)'^ji^ The pulses are displayed in the top of Fig. 1, with the different width 
A® = 4®,„,4/3®,„,0.8®,„,0.4®„,. 

One has the increased entanglement following the increase of pulse intensity. The entangle¬ 
ment keeps increasing though the pulse intensity has started to decrease. Only until the pulse 
almost disappears, does the entanglement begin to decrease and oscillate for a period. The en¬ 
tanglement will exist longer than the pulse duration. 

The entanglement under blue detuned central frequency [Fig. 1(a) and Fig. 1(b)] is obviously 
larger than that under red detuned central frequency [Fig. 1(c) and Fig. 1(d)]. However, a main 
difference of the pulse driven entanglement from the CW counterpart is that a spectrum of fre¬ 
quencies with disparate detunings give rise to the actual entanglement. As seen in Fig. 1, the 
entanglement values for the pulses of narrower bandwidth behave monotonously in both blue 
and red detuned regime for the central frequency, while those of larger bandwidth show the 
entanglement patterns of both blue and red detuned frequency components. The entanglement 
from red detuned frequency components oscillates with time li25]l . so this pattern always ex¬ 
ists for the narrowest pulse (corresponding to the widest frequency bandwidth) in Fig. 1. As 
the result, the entanglement for the pulses of larger bandwidth show insignificant difference 
in Figs. 1(a) and 1(b), except that more oscillation pattern is introduced in 1(b) as the central 
frequency moves toward the red detuned regime. When the pulses’ central frequency becomes 
red detuned, the oscillation for one of the pulses will disappear due to the suppression of en¬ 
tanglement for its whole frequency spectrum in the regime. On the other hand, given the same 
driving intensity, the entanglement will become higher with a narrower width. 

For comparison, we consider other pulses of different profiles, such as square pulse, triangle 
pulse, sawtooth pulse and trapezium pulse; see the top of Fig. 2. Their evolutions for the cen¬ 
tral frequency detuning Aq/®,,, = — 1 [Fig. 2 (a)] are similar, so the different pulse shapes do 
not affect the entanglement much in this case. Meanwhile, difference appears in the evolution 
with red detuning Aq/®,,, = 1 [Fig. 2 (b)]. Compared with Gaussian pulse driven OMS, the 
entanglement values under the pulse drives in Fig. 2 are higher around the SQ resonant point 
(Ao/®,„ = —1). It can be explained by the fact that a Gaussian pulse contains more significant 
off-peak frequency modes, and the frequency components whose detunings are far away from 
SQ resonant point contribute less to the entanglement. This is also the reason for why a nar¬ 
rower bandwidth of Gaussian pulse is better for entanglement creation. Moreover, due to their 
more centered frequency spectra, the differences of the entanglement evolutions in the blue and 
red detuned regime become more obvious for the pulses in Fig. 2. 

5.2. Entanglement evolution under a series of laser pulses 

As show in Fig. 1, the entanglement could last for longer time than the pulse duration. The 
entanglement may survive if one more pulse drive the OMS before it disappears. Since the 
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Fig. 1. Entanglement evolution under Gaussian pulse drive with different widths (displayed 
on the top). Long dashed (green) curve is for Aco/(b„, = 0.4, short dashed (black) curve 
for Aa)/co„, = 0.8, thin solid (red) curve for Atu/o),,, = 4/3, and thick solid (blue) curve 
for Ato/co,,, = 4, respectively. The corresponding detuning Aq/co,,, for Figs. 1(a)-1(d) are 
A()/cOi„ = —1, —0.5,0.5,1, respectively. The other parameters aieg/ K = 10^®, (Oj„/k = 2.5, 
f^mhm = 10^, and T =0. 
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Fig. 2. Entanglement evolution under pulse drive with different profiles (displayed on the 
top), (a) The detuning is Aq/co,„ = — 1. (b) The detuning is Aq/ to,,, = 1. The other parame¬ 
ters are the same as in Fig. 1. 


entanglement will be kept and raised in this case, the minimal value of entanglement will be 
nonzero and the peak value of entanglement can be slight larger than that under a single pulse. 
This phenomenon will occur if the gap time between the centers of the pulses is not too large. In 
Fig. 3, we show the entanglement evolution of the OMS driven by a series of Gaussian pulses. 
The gap time (A(Kt) = 8) is relatively large in Fig. 3(a), resulting in the repeated evolution 
pattern of the entanglement under single pulse driven OMS [see Fig. 3(b)]. When the gap time 
(A{Kt) =5) is relatively small as in Fig. 3(c), the entanglement evolution will be totally different 
as shown in Fig. 3(d). The peak value is about 1.6 for the width Aw/®,,, = 4 and 1.1 for the 
width Aco/cOm = 4/3, which are slightly larger than the ones (1.5 and 1, respectively) shown 
in Fig. 1. At the same time, the minimal value becomes nonzero (about 0.1) for the width 
A®/®,„ = 4. Especially, for the width A®/®,„ = 4/3 (the overlap of two pulses can be neglected 
in this case), the minimal value will be even larger (about 0.3). A rather stable entanglement 
(> 0.3) can therefore be achieved by pulsed drive, even the overlap of two pulses is not so 
obvious. 

5.3. Entanglement evolution under high temperature 

In the above discussions, the temperature is set to be T = 0. However, temperature is also a main 
factor that affects the existence of entanglement. In Fig. 4 we show the entanglement evolution 
under the condition = lO'^ while other conditions are the same as in Fig. 1. Obviously, the 
entanglement decreases with high temperature, and the phenomenon of entanglement sudden 
death and sudden revival becomes apparent. Through the comparison, one sees that the robust¬ 
ness of the entanglement against high temperature for the blue detuned frequency components 
















Fig. 3. Entanglement evolution under a series of Gaussian pulses with the detuning 
Aq/co,,, = — 1. The left two figures are for the pulses and the right ones are for the en¬ 
tanglement evolution. The gap time of the pulse centers in (a) and (b) are A{Kt) = 8 and 
A( Kt) = 5, respectively. The width for thin (blue) curve is Am/®,,, = 4, while that for thick 
(red) curve is Am/m,,, = 4/3. The other parameters are the same as in Fig. 1. 


is much more significant than that from red detuned ones. As it is also shown in these figures, 
the entanglement evolutions under higher temperature become more complicated, because they 
result from more intricate interplay of various physical factors (temperature, frequency band¬ 
width, central frequency detuning and drive intensity). The competition between the factors 
gives rise to the seemingly unexpected entanglement amplitudes in Figs. 4(b), 4(c) and 4(d). 
At the central frequency detunings in 4(b) and 4(c), the higher driving intensity of the pulse 
with the widest bandwidth actually suppresses the entanglement value, because the thermal 
noise effect killing the entanglement is enhanced by the driving strength [see Eq. (15)]. If the 
central frequency is exactly BS resonant as in 4(d), the suppression of entanglement for the 
whole spectrum becomes more significant with the temperature, leaving the remnant due to the 
far-away frequency components from the point, which is only possible for the pulse with the 
widest bandwidth. To achieve high and robust entanglement, blue detuning around Aq/co,,, = — 1 
associated with narrow pulse bandwidth, should be the best choice. 

6. Conclusion 

We have depicted the optomechanical entanglement evolutions under various pulsed drives. 
Compared with CW drives, high and robust entanglement can also be created especially with 
pulses of blue detuning and narrow bandwidth. Such entanglement can last for long time with 
repeated series of driving pulses. The flexibility in engineering pulse frequencies and shapes 
enables one to achieve variety of entanglement evolution patterns for OMSs. 
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Fig. 4. Entanglement evolution under Gaussian pulse drive with the temperature n,„ = 10^. 
The other parameters are the same to Fig. 1. 
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